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Internetting of Fires 


by David Olwell and Alan Washburn 


Abstract 


The advent of long-range, accurate weapons, together with advanced sensors and information systems, 
has changed the nature of the modern battlefield. Targets can now be engaged by weapon systems that are 
widely separated geographically from their targets. These new capabilities bring increased demands for 
coordination and optimization when multiple weapon systems (“shooters”) engage multiple targets. There 
are aspects of an assignment problem where each shooter does what he is good at relative to other shooters 
in order to maximize destruction to any given set of targets. This report introduces methods for making 
assignments optimally, and compares them to simpler methods that optimize only locally, instead of 
globally. The benefits for global optimization depend on circumstances, but can be substantial. 


Executive Summary 


Here we introduce a simplex-type algorithm for optimally assigning shooters to targets, and 
compare it to other methods. The algorithm rapidly solves the problem of optimally assigning a 
fixed collection of shooters to a fixed collection of targets, given known kill probabilities for 
each potential assignment. The algorithm itself is represented by SJMPLL.dll, a Windows 
dynamic-linked library. 

Optimal methods, of which SIMPLL is but one example, are compared to other methods that 
make decisions locally, rather than globally. These comparisons are made using NetFire.xls, an 
Excel workbook that implements all methods simultaneously so that they can be compared on 
specific eenanos Optimal methods are bound to win such a comparison, but the interesting 
question 1s whether the margin of victory is big enough to justify the C4ISR system that must 
accompany the global point of view in practice. The answer, unfortunately, is that the margin of 
victory depends strongly on circumstances. In cases where all targets are equal and all shooters 


have the same capability, there is little to be gained by global coordination. In more variable 


circumstances, the margin of victory can be large. 

It is sometimes suggested that warfare can be conducted by autonomous “agents” that 
operate and make decisions based only on locally available information. It is also sometimes 
suggested that modern armed forces should attempt to achieve a C4ISR system sufficiently 
sophisticated to achieve “internetting-of-fires”, since it is only with such a system that the 
capabilities of long-range weapons can be exploited to the fullest. These suggestions are in 
conflict, since the one asserts that local optimization is sufficient or at least inevitable, while the 
other asserts that only global optimization will do. In an abstract way, these ideas are 
represented by the different targeting schemes discussed in this report. Possibly NetFire.xls 
and/or SIMPLL.dlil will prove useful to those interested in exploring this issue further. 

We have found instances where global use of information can achieve the same results as 
local use of information, but at a cost of only half the number of weapons. We have also found 
instances where there is no advantage to global use of information. We explain some of the 


factors that influence the difference. 


1. Introduction 


Modern warfare is distinguished by the presence of increasing amounts of relevant 
information on the battlefield. The fog of war will never disappear, but it has thinned a bit due to 
continuing improvements in sensors, computers and communication networks. It 1s now realistic 
to talk about “internetting of fires” in the sense that targeting information used by a given shooter 
may come from sensors not under his direct control. But this internetting capability comes at a 


stiff price. There are even questions about whether the money might be better spent in more 


effective ways. General Howell M. Estes II], USSPACECOM, has put it succinctly: 











Hard choices need to be made between investments in information infrastructure [and] 
the combat systems themselves. This is an extreme dilemma, because combat systems, 
without timely, relevant information, are useless. On the other hand, you cant take out 
an enemy tank with just information. We need to strike a balance between ‘shooters’ and 
‘information systems’ if we’re going to be successful in the future. 


The subject of this report is the value of information in the context of land combat, with 
emphasis on the tradeoff between information and firepower. To get at this tricky issue, we 
postulate a group of targets, each of which has a value, and compare various targeting schemes 
on the basis of the Measure of Effectiveness (MOE) “fraction of total value killed, on the 
average”. Clever schemes with a high investment in information systems can be expected to 
achieve high scores, but so can less advanced schemes that simply have more firepower. The 
incremental dollar should be spent in whatever manner increases the MOE the most. The 
resulting combat model, NetFire, is embodied in an Excel workbook. 

The general concept in NetFire is that there is a certain length of time called a decision cycle 
within which no Battle Damage Assessment (BDA) is possible. Within that cycle, a fixed set of 
weapons (“shooters”) is available to be fired at the enemy. NetFire reports results for several 
different methods of targeting, some methods requiring only local information and some 
requiring global information. One of the latter is a system wherein each shooter is directed to 
fire in such a manner as to maximize the MOE. The algorithm for doing this (SIMPLL) 1s a 
principal output of this research, and may be of use more generally than in NetFire. 

NetFire is an abstract combat model. The enemy does not move or shoot back, and, with one 
exception, time is not represented. “Combat” simply consists of one-off situations where 
dispersed shooters must simultaneously commit to targets. Nonetheless, NetFire can be used to 
explore the significance of certain kinds of information, specifically BDA information and 
coordinating information. The next section describes the analysis that underlies the calculations. 


The specifics of using the Excel workbook NetFire.xls will be described later. 





2. Analysis of Fire Allocation Methods 
Cases RANDOM, GREEDY, and OPT 

Information may or may not be important to determining the battle outcome. Consider the 
following scenario. Suppose that there are m shooters and n targets, with each shooter having a 
single shot that is capable of killing its target with probability Pk. Temporarily make the 
following assumptions: 

e the weapons are all “long range” in the sense that any shooter can physically attack 
any target; 
e each target has the same value; 
e each shooter possesses the required targeting information for each target; and 
e there is no BDA between shots over the time period under consideration, so all shots 
might as well be committed in one salvo. 
It does not follow from these assumptions that shooters can communicate with each other, nor 
does it follow that there is a commander capable of sensing the overall situation and giving 
appropriate orders. How much is this additional coordinating information worth? 

The worst possible outcome from the viewpoint of the shooters would be if all m shots were 
directed against a single target, in which case the unfortunate target would be overkilled (unless 
Pk is very small), while n-1 targets would not be attacked at all. The optimal outcome in this 
case would be if all targets were shot at equally. Clearly, achieving the best outcome requires a 
commander with perfect information. However, the worst outcome a/so requires a commander 
with perfect information, albeit with a different motivation, so comparing the best with the worst 
will not be enlightening. If the shooters cannot communicate—either with each other or with the 


commander—then it makes sense to assume that the target chosen by one shooter is statistically 








independent of the targets chosen by other shooters. The resulting outcome will be intermediate 
in terms of outcome, but represents an information extreme in the sense that no coordination 
whatever is required to achieve it. It is these two extreme situations, random — optimal, that 
should be compared in discussing the value of coordinating information. 

The optimal allocation will allocate the shots equally, so m/n shots will be assigned to each 
target. Each target survives only if all of these shots miss, so the average number of targets 
killed will be £,,, =n(i-(- Pk)”"”). On the other hand, if each shot is directed at a random 
target, then each target survives a given shot with probability 1-Pk/n, independently of the other 
shots. Since there are m shots in total, the average number of targets killed will be 


E 


rand 


=n(1-(1—Pk/n)”). Figure 1 shows both of these quantities as a function of m when n=10 


and Pk=0.5. 
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Figure 1: Average number killed versus total number of shots. 


It can be seen from Figure 1 that there is very little difference between the optimal and 
random results when m is zero (in which case both quantities are of course 0) or very large (in 
which case both approach n). Even in the middle, where mPk and n are approximately equal, the 


functions are surprisingly close, given the extreme differences between the command and control 


systems being modeled. When all weapons are equally effective, all targets are equally valuable, 
and the kill probability is not too large, there turns out to be only a limited payoff for information 
that would permit the commander to make an optimal assignment for every weapon. When the 
global object is simply to distribute weapons evenly over targets, the local tactic of aca saing 
weapons uniformly at random is not that far from global optimality. 

In mixed scenarios where target values differ and the kill probability depends on both 
weapon and target, the situation changes. It is no longer clear what firing strategy should be 
used 1n the extreme case of no coordination, since now there is a conflict between treating all 
targets equally and selecting the target for which the kill probability is largest. To be more 
precise about this, as well as to account for the type of internetting information that makes firing 
even feasible, we introduce the following notation: 

i indexes shooters, & indexes targets 

V,, = value of target k 

S==number of shots from shooter i 

Dix = prob that one shot of 7 kills k, given i shoots at k 

I= 1 if it is feasible for i to shoot at k, else 0 {depends on weapon range, internetting, and 
possibly other command and control factors} 

Pi=pitlix {only the product of p and J is important} 

When there is no commander able to coordinate fires, each shooter could still choose a target 
at random. But a target truly chosen at random could very well have a zero kill probability, so 
instead we assume that the randomization is done only over those targets for which the kill 


probability is nonzero, the ones that are reasonably close to the shooter. To be precise, we define 


Case RANDOM: Assume that each shot from 7 is equally likely to be taken against any 











target k for which P;>0. Define the sets T={k|P;>0} and R= {i|/P;>0}, and let #() count the 


number of elements in a set. Then P;;/#(T;) is the probability that each of 7’s independent shots 


kills target k, and therefore the probability that target & survives is gq, = I] ( _ A, | . The 
iéR, i 


survival events are not independent among the targets, but nonetheless the expected remaining 


target value is V*= > Vi(i-4q,). 
k 


Case RANDOM could still be criticized for sometimes directing fire against targets where 
the kill probability is low. It is tempting instead to assume that shooter i directs all of his fire 
against the target for which the product of target value and kill probability is largest. This 
“GREEDY” procedure is also feasible without coordination by a commander as long as target 
values can be sensed locally. To be precise, we define 

Case GREEDY: Assume instead that each of i’s shots is taken against the target for which 


it is most effective, namely target K(i)=argmax;,{V;Pix}, and let R= {ilK(i)=k} be the set of 


shooters who attack target A. Then q, = [Ia — P.)', and the expression for * is as before. 


ieR, 

Case GREEDY may not do as well as Case RANDOM. Suppose that there were some single 
valuable target that is very vulnerable to fire from every shooter. Then all fire would be directed 
against that one target in Case GREEDY, and all others would escape. Case GREEDY makes no 
provision for spreading fire over other targets. 

The fact that neither RANDOM nor GREEDY is clearly the right thing to do when no 
coordination is possible simply bespeaks a difficult fire allocation problem. In this circumstance, 
one might suspect that a coordinating commander could make more of a difference, particularly 


when at least some of the kill probabilities are large. Our reasons for introducing the RANDOM 





and GREEDY procedures are not that we advocate them in any sense, but that they are simply 
two procedures that might characterize what actually happens in the absence of the coordination 
that internetting makes possible. Specifically, they might characterize target selection by 
“agents” whose behavior is based only on locally sensed conditions. 

The best scenario for coordination information, then—the case where its value in terms of 

additional targets killed is largest—will be when: 
e targets have different values; 
e kill probability is variable; 
e at least some kill probabilities are large; and 
e the scenario 1s not one of massive overkill or underkill. 

It might seem that this is a rather narrow scenario, and in a way it is. Nonetheless, our 
intention is to concentrate on cases of this sort. Our contention 1s that this situation characterizes 
modern warfare. In particular, modern weapons may be guided weapons that are long range and 
have large kill probabilities, but which are in limited supply on account of being expensive. It is 
particularly important to employ such weapons carefully, reserving them for cases where less 
_ expensive weapons cannot do the job. Only coordinating information can enable this. 

When coordinating information 1s available and timely commands to individual shooters are 
feasible, the fires of all the shooters can be directed in a manner that maximizes the average 
target value killed. The resulting optimization problem might be called a nonlinear assignment 
problem, since the objective will be best accomplished by assigning each shot to the target at 
which it is most effective. Let the variable x; represent the number of shots made by shooter i at 


target k. Then Case OPT 1s to 


maximize DY, (1—q,), subject to the constraints 
k 











qq, = | ]a-4.* ) 
DG <S,, 

k 
Xx, 29 
Here, and in all following optimization formulations, we use the convention that variable names 
begin with lower case letters, whereas constants begin with upper case letters. The formula for 
qx, the miss probability for target k, reflects the usual independence assumption in the extended 
product of miss probabilities. 

Case OPT is actually a nontrivial Mathematical Programming problem, since the objective 
function is nonlinear and the variables are also required to be integral in value. We cannot solve 
it in reasonable amounts of time for large problems. However, an upper bound on the optimal 
objective function can be obtained by relaxing the integrality constraint, which produces a 
simpler, continuous optimization problem for which an efficient solution technique is known 
(Washburn[1995]). A lower bound can then be constructed by rounding the solution to integers 
in such a manner that shot constraints are still enforced. This process 1s carried out in 
SIMPLL.dil, a dynamic-linked library that can be used with any Windows-based application. 
The SIMPLL algorithm is actually a branch-and-bound method (Washburn[1998]) that will find 
a solution accurate to within TOL, a nonnegative tolerance for error. 


Other Possible Views of the Fire Allocation Problem 


The basic viewpoint taken above is that finite stocks of weapons must be allocated to targets 
within any given decision cycle. The weapons themselves are not “costly” in any sense other 
than being in limited supply. This is a common point of view, but by no means the only 
reasonable one. If weapons are costly, instead of being constrained, it may make sense to 


minimize the cost of “getting the job done” in some sense. This new point of view results in the 


10 





consideration of optimization problems such as NLP1: 


minimize Cx , subject to the constraints 
i,k 


a =[]G-7)*, 

> Vg, SV, 
k 

X, 20 

In NLP1, the weapon constraints S$; have been replaced by weapon costs C;, and the object is 
to minimize the total cost of shooting, subject to the total surviving value not exceeding V. Costs 
may be in units of dollars, or may reflect the weight or volume of the weapons, which would be a 
consideration when moving weapons into a distant theater. Cost might also reflect tactical 
preferences. 

NLP1 is a nonlinear but convex program, since the target value constraint amounts to 
requiring that a convex function not exceed V. NLP1 can therefore be solved by any of a number 
of commercially available nonlinear optimization packages. If the variables are required to be 
integers, a branch-and-bound algorithm could be based on an NLP1 relaxation, as in SIMPLL. 
However, there is no tailored algorithm for solving NLP1. We have investigated solutions to this 
problem using GAMS (Brooks, Kendrick, and Meeraus[1988]) to access general purpose solvers, 
and find that it solves much more slowly than Case OPT. For example, solving a 25 shooter, 75 
target problem takes about a second for SIMPLL, but a minute for GAMS using the formulation 
shown in the Appendix. The performance spread increases as the problem size grows, 
highlighting the advantages of having a special-purpose tailored algorithm such as SIMPLL. 


In principle, one could also solve NLP2 


minimize > Cz; +> V4, , subject to the constraints 
i,k k 











qk = [Ia ia 

Xe = 0 
NLP? is the problem of minimizing total cost, including the “cost” of the surviving targets, 
subject to no constraints on either target value or weapon consumption. Weapon consumption is 
left to be whatever is natural for the target complex. NLP2 is a Lagrangian relaxation of Case 
OPT as well as of NLP1, so it should be easier to solve than either one. A special-purpose 
algorithm similar to SIMPLL could surely be developed for it. NLP2 suffers from the practical 
problem that both weapons and targets must be valued on the same scale, probably dollars. The 
tradeoffs that are implicit in Case OPT and NLP1 are explicit in NLP2. 

One could also go in the opposite direction and add constraints to NLP1, rather than remove 
them. If every target is constrained to survive with only an input probability Q;, one arrives at 
NLP3: 


minimize C,x,, , Subject to the constraints 
i,k 


q, = [la = are ’ 


q, SQ, 
x, 20 


NLP3 is especially simple because q; can be eliminated, after which the constraint set is: 

> *x nd — F,)<In(Q,) 

420 
Since the objective function and all constraints are linear, this is a Linear Program. The problem 
here is that one must give a survival probability for each target, with no tradeoffs possible 
between targets. The minimized objective function might be disappointingly large if QO; were 


made small for some hard-to-kill target. 
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Case BDA 


It was mentioned in the Introduction that NetFire does not model time, with one exception. 
Here we discuss that exception. 

All of the previous calculations pertain to a “one-off” war that is one decision cycle in length, 
where all weapons are assigned without regard to the effects of other weapons. It is likely in 
such circumstances that some weapons will be assigned to targets that have already been killed 
by other weapons—a potentially significant waste if the weapons are expensive. In a multi-cycle 
war where BDA information about the status of targets is available, this waste could be 
prevented. To bound the value of such BDA information, we consider in this section a situation 
where there is no time pressure, and where every shot is followed by an instantaneous, perfect 
assessment of whether the target was killed or not. No real BDA system is this good and no war 
is so leisurely, but it is still of interest to see how much additional target value could be killed (by 
weapons not previously wasted on dead targets) in this circumstance. 

The data Pix, S;, and V; has the same meaning as before, but many policies are now feasible 
that were previously not. These new policies all involve shooting a single shot at some target, 
assessing its effects, and then, depending on the effects, shooting another single shot at some 
possibly different target, and so on. Each of these policies is a very complicated object, since all 
eventualities sate be allowed for. Each might be expressed as a decision tree. For example, one 
policy might begin, “Fire weapon 3 at target 6. Then follow branch A if target 6 is killed, 
otherwise follow branch B”. There are a great many possible such policies, but the number is 
finite. The problem of finding the one that maximizes the expected target value killed is 
therefore well defined, but the number of possible decision trees is large enough to make 
exhaustion an unattractive technique. 


We know of no practical way of finding the optimal policy. Our goal here is merely to 
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bound from above the expected target value that such a policy might achieve. Toward that end, 
we define a collection of indicator random variables: 

Y;, =] if target k is killed by weapon i (for each k, at most one of these can be 1), 

Y,=1 if target k is killed, and 


Xit=1 1f weapon i is assigned to target k (for each 7, at most one of these can be 1). 


i 


The firing policy will induce many correlations between these random variables, so 


independence assumptions among them are not appropriate, but the collection is still useful for 


formulating the problem of finding the optimal policy. We first note that Y, = Pe , and that the 


total value killed is Z = YW, . The problem (call it P1) of finding the optimal policy can 
k 


therefore be posed as maximizing z, the expected value of Z, subject to constraints 


(1) >) X, <5, for all i with certainty, 
k 
(2) Y,= pee for all & with certainty, 


(3) Y;<1 for all & with certainty, 


(4) other constraints. 


The other constraints in P1 include the crucial relationship between Xj, and Yj, the essence of 
the firing policy. Whatever that relationship, it must certainly hold that E(Yin)s Pix E(Xx ). Now, 
using lower case letters for expected values of random variables, we can construct a relaxation of 
P1 that we name P2: 


maximize z= y Vy, Subject to 
kK 


(1) ye < S, for all 7, 
k 





(2) », S >) x,P, for all k, and 


(3) y,< 1 for all k. 

P2 is a relaxation of P1, because a relationship that is true with certainty will also be true on 
the average, because sums and expected values can always be interchanged, because (4) implies 
that E(Vx)<%ix Pix , and because the rest of (4) has simply been omitted in P2. P2 is a simple 
Linear Program that provides an upper bound on what is achievable with any shoot-look-shoot 
policy. P2 has a direct interpretation where x; is the average number of weapons of type i used 
on target k, and y,; is the probability that target k is killed. (1) requires that the average number of 
weapons of type i used not exceed the number available, (2) requires that the effect of each shot 
not exceed P;,, and (3) requires that target k be killed at most once. 

This BDA upper bound is one of the computations made in NetFire.x/s, using Excel’s Solver 


to solve the Linear Program. 


The AO Boundary 


Hierarchic military organizations sometimes partition the target set into parts, with each unit 
in the hierarchy being responsible for one of the parts. Doing so simplifies command 
relationships and minimizes the need for communications. It may also be beneficial in avoiding 
some of the overkill errors that Case GREEDY is subject to, where too many shooters choose a 
soft or valuable target. However, the partition can also be viewed as the addition of constraints 
to Case OPT, which can only have the effect of reducing the objective function. Is the reduction 
significant? An easy way to find out is to partition the battlefield into two halves, with the 
shooters and sensors in one half being ineffective in the other. Each of the cases discussed above 
can be applied separately in each half, and the results added. By comparing the results with the 


cases where the battlefield is not partitioned, one can measure the significance of the partition, as 








in NetFire.xls. 


3. NetFire.xls 


JANUS Data 


NetFire is distributed with a variety of defaults, including the entire database pertaining to 
(shooter, target) kill probabilities. This subsection will be of interest only if the user wishes to 
change those defaults, and may be skipped on first reading. 

Most inputs to NetFire.x/s are directly from the user. The main exception is a mechanism for 
incorporating kill probability data from the JANUS combat model. To do so, one must first 
obtain from the JANUS database a text file named PX/L.csv in the same format as sheet “PKIL”. 
There should be a single header line followed by a separate line for each (weapon type, target 
type) pair, typically a long file because each shooter type may have several target types in 
JANUS. Each line after the first should give a sequence of kill probabilities (more precisely “hit 
and kill” probabilities) indexed by range and attitude relative to the shooter, as shown on sheet 
“PKIL”. This file should be placed in the same directory as NetFire.xls, after which pressing the 
“Update tables from PKIL.csv” button on sheet “PKIL” will overwrite the data on that sheet with 
the new data in PKIL. csv. 

The command button labeled “Update tables from this sheet” on sheet “PKIL” uses the data 
on that sheet to create sheets “QTable” and “RTable”. As NetFire is distributed, the data shown 
on sheet “PKIL” is abstracted in subroutine NewQK( so that every (shooter, target) pair has a 
single miss probability that is valid out to a maximum range, after which the miss probability is 
1. This subroutine then writes the miss probabilities to sheet QTable and the maximum ranges to 
sheet RTable, each in the form of a table with shooters on the rows and targets on the columns. 


Making miss probability be (say) a linearly interpolated function of range would require writing 
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new VBA code in subroutine NewQK(). 

Subroutine NewQK() also incorporates a weapon reliability parameter that is hardwired to .9. 
As aresult, no miss probability on sheet QTable will be smaller than .1. This subroutine also 
initializes the number of shots available (the default for column 1 is 4 shots) and the sampling 
weight (the default for column 2 is set to 1 for all shooter types) of each shooter, as well as the 
values (row 1 is set to 10) and sampling weight (row 2 1s set to 1) of each target. These values 
can subsequently be edited by hand, but will be renewed at the default values whenever the 
“Update tables from this sheet” button is pressed. The defaults can be changed by editing the 
VBA code for NewQK(). 

As distributed, NetFire.x/s is unclassified because it uses an unclassified dataset from 
JANUS. More realistic, classified data would provide more realistic, classified results. 


The Main Sheet of NetFire.xis 


The main sheet of NetFire.x/ls is the sheet named “pij”. On this sheet can be found a digital 
and graphical description of the current scenario, along with command buttons “Resample’’, 
“ValKill’, and “MultiRep”. These buttons have the following functions: 

e The battle involves three types of objects: shooters, sensors and targets. Command 
“Resample” randomly selects and writes the object properties to the “pi” sheet, 
thereby defining a scenario. It also shows the object locations in a chart. The user’s 
inputs here are the cells colored with a light green background: the number of each of 
the three types of object, the dimensions of the areas within which they are located, 
and, if desired, an AO boundary location. Resample also uses data from the sheets 
“QTable” (miss probabilities), “RTable” (max ranges), and “STable”’ (sensor data) to 
compute and write the miss probability matrices for both the AO case and the non- 


AO case. The user input cell “Value variance” increases the variability of the target 


ee 











values, with input 0 leaving them all at the default. 

e “ValKill” computes the MOE “average fraction of target value killed” for four cases 
and two AO possibilities, writing values to the rows and columns of the output cells 
with a light yellow background. One method (OPT) has two rows because both a 
feasible solution and an upper bound need to be shown. ValKill uses the data in the 
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ranges named “shots”, “values”, “missprob”, and “missprobr!”. Two tables of miss 
probabilities are needed because results must be computed both with and without an 
AO. This data can be edited directly if desired, but will be overwritten when the 
Resample button is again pressed. 

e ‘“MultiRep” iterates the Resample/ValKill combination an input number of times, 
leaving only the last replication on sheet “pij”. The 10 numbers computed by ValKill 
are written on the sheet “MultiRep” in rows, plus one additional row that averages the 
results in each column. The desired number of iterations is in the green cell right 
above the button. 

Figure 2 shows a screen shot of the main sheet, including part of the map showing object 
locations and two of the command buttons. The check boxes are included because sometimes it 
is useful to resample only parts of the scenario, but normally all will be checked. The 10 
numbers shown next to the ValKill button are the main outputs. 


Scenario Development in Resample 


The Resample command button creates each scenario by first creating objects called 
shooters, sensors and targets with randomly selected properties. All such objects have locations 
(.X and .Y) that are uniformly sampled based on the limits given by the user on sheet “pi”. 


These coordinates are shown on that sheet both digitally and graphically. 
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Figure 2: Screen shot of the main sheet “pij” of workbook NetFire.xls. 


Each shooter also has a type that is selected with a sampling weight read from column 2 of 


sheet “QTable”. These weights are all set to 1 when the sheet is created, but can be changed by 


the user if some shooters are desired to be more numerous than others. For example, make the 


weight of “PISTOL” 10 if Pistols should be 10 times as numerous as other shooters (only relative 


values of weights are significant, so arbitrary nonnegative values are permissible). Other 


editable items taken from “QTable”’ include: 


e The number of shots available to each shooter (column 1). 
e The value of each target (row 1). 


e The sampling weight of each target (row 2). 
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e The miss probability for each (shooter, target) pair, conditional on a shot being 
feasible. 
The maximum range for any given (shooter, target) pair is taken from sheet “RTable”. 

The Resample command also creates sensors, using data on sheet “STable”. Each sensor has 
a name, a sampling weight, and a range within which it detects and identifies all targets. Given 
the user input number of sensors, the sensors are chosen randomly according to their weight, 
assigned the appropriate name and maximum range, and located uniformly at random within 
user-designated limits. As with shooter and target objects, the properties of each selected sensor 
are shown on sheet “‘piy”’. 

Finally, the Resample command creates two miss probability tables and shows them on sheet 
“pi”. The first table is located in named range “missprob’’, and applies to the case where all 
targets are in the same area of operations (AO). The miss probability for any individual (shooter, 
target) pair 1s 1 if the shot is not feasible, or otherwise the value read from sheet “QTable”. 
Feasibility requires that: 

e The target must be within firing range of the shooter. 
e The target must be known. A target is known if it is either 
o within range of any sensor, or 
o within VisRange of the shooter. 
The VisRange parameter is hardwired at 3,000 meters in subroutine Sample(), which is the 
subroutine called by the Resample command to do all of the scenario generation work. 

A separate table “‘missprobrl” is created for the AO case. All objects are partitioned by the 

AO boundary, and no shooter on one side of the boundary can even learn of a target on the other 


side, much less shoot at it. The AO miss probabilities are always larger than the non-AO miss 


20 





probabilities, since shooters and sensors in one AO are not allowed to report on or shoot at 
targets in the other. It does not follow that results in the AO case will always be worse (see 
discussion below), although they usually are. 


Computations Carried Out Within ValKill 


All of the inputs required by the ValKill command are displayed on sheet “‘pij’”’. In fact, the 
only data involved are in the named ranges “shots”, “values”, and “missprob” (non-AQ) case, 
and “missprobrl” (AO case). All of this data can be edited directly before pressing the command 
button, but it will be overwritten the next time the Resample command button is pressed. 

The function of the ValKill command is to compute the 10 numbers in the yellow shaded 
area, all of which are the MOE “average fraction of target value killed” in different 
circumstances. All of these computations are carried out within the ValKill() subroutine. The 
two local firing rules, RANDOM and GREEDY, are implemented directly within the subroutine. 
The OPT firing rule is implemented by calling subroutine SIMPL() in the external dynamic- 
linked library SIMPLL.dll. This subroutine returns an upper bound on what is achievable by any 
firing procedure faced with the same problem, so the upper bound 1s displayed, as well as the 
score achieved by the (nearly) optimal firing plan. The subroutine is declared at the beginning of 
module 1 of VBA code, and called from subroutine ValKill(). Finally, the BDA upper bound 1s 
calculated by calling subroutine LP _Upr(), which formulates and solves the appropriate Linear 
Program using the Excel Solver and sheet “LP Upper’. Solver will only handle problems with 
200 or fewer variables. If the product of the number of shooters and the number of targets 
exceeds this limit, the upper bound is simply reported to be zero. This constraint could be 
eliminated by replacing the standard Solver with a better optimizer. 

Each of the computations is repeated in the AO case, except that miss probabilities are taken 


from “‘missprobr!” instead of “missprob”’. 
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The SIMPL()subroutine in S7MPLL.dll can theoretically be called from any Windows 
application, as long as declarations are made carefully and the Compaq FORTRAN runtime 
library DFORRT.dll 1s also available. 

NetFire was designed to generate random scenarios for analysis. It is also possible to use 
NetFire to analyze existing scenarios. To do so, first adjust the numbers of targets, shooters, and 
sensors to the correct values and press the Resample button. This will generate a random 
scenario with the correct dimensions. Then enter the correct target, shooter, and sensor data 
manually into the appropriate rows and columns. Save the workbook to protect all this work. 
Then uncheck the boxes for targets, shooters, and sensors under the Resample button on page 
“pij”. Pressing the Resample button will then calculate the two required tables of miss 
probabilities. The ValKill button will then evaluate the different firing strategies. 


General Tendencies 


There are a few theoretical inequalities that must hold regardless of scenario. These will be 
stated in terms of VK(i,), the i” row and the a column of the 10 outputs, i=1,..,5; j=1,2. Table 1 
shows a typical example. 

1. VK(i,j)S VK(4,j); i=1,2,3 and j=1,2. This is because an upper bound on the optimal 
score must necessarily be greater than any feasible score, and all scores shown in the 
first three rows are feasible. 

2. VK(i,2)< VK(4,1); i=1,2,3. VK(4,/) is an upper bound on what can be achieved by 
any firing method that is required to make all shots within one cycle when there 1s no 
AO boundary, and therefore is also an upper bound when the AO constraint is added. 

3. VK(4,)SVK(5,/); j=1,2. One basically gets to the BDA bound by replacing the 


function 1-exp(-x) in the OPT bound by the dominant function min(1,x), so the BDA 


ZZ 


bound is necessarily at least as large as the OPT bound. 





Single AO Two AOs 
ValKillRand 0.2552852 0.245430418 
ValKillGreed 0.26754648 0.245757911 
ValKillOptLwr 0.31123148 0.276201159 


ValKillOptUpr 0.31165106 0.276259909 
ValKillSLSUpr 0.3365397 0.29697468 


Table 1: Typical outputs from NetFire.xls. 

























All of the above inequalities should be true regardless of scenario. If any should ever fail to 
hold, then please notify one of the authors because a bug has been revealed. 

Other than these inequalities, not much can be said in general about even the ordering of 
results, much less about absolute differences between them. It is not always true that GREEDY 
is better than RANDOM. To provide a counterexample, generate a scenario where all miss 
probabilities are equal. Every shooter will fire at the most valuable target in the GREEDY 
method, whereas the RANDOM method will spread fire equally over all. If the most valuable 
target is only slightly more valuable, RANDOM will get a better score. Also, in spite of 
inequality 2 above, it is not always bad to have two AOs. The GREEDY method, in particular, 
may do better when there is a constraint that prevents certain shooters from shooting at certain 
targets. 

The incremental value of information is highly circumstance dependent, and there are so 
many dimensions to a NetFire scenario that one hesitates to make general conclusions beyond 
the above inequalities. Nonetheless, after generating and solving a large number of scenarios, 
any user will come to conclusions about tendencies. For the scenarios generated by the authors 
in the course of developing NetFire, the following tendencies have been noted: 

e GREEDY generally outperforms RANDOM, but not by much. 


e The OPT policy typically increases the fraction of target value killed by about 20%, 
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relative to either GREEDY or RANDOM. 
e The BDA upper bound 1s typically about 30% above the GREEDY/RANDOM level. 
e A few percentage points are usually lost in going from one to two AQOs. 

The performance of BDA is surprisingly low, since BDA as used here represents a perfect 
assessment system where time is not constrained and the true status of every target is always 
reported correctly after every ie As a means for promoting efficiency 1n Aung: BDA appears 
to be not particularly effective, at least 1f the scenarios examined by the authors are typical. Of 
course, knowing the correct status of targets has a military payoff beyond weaponeering 
efficiency. 

All of the above observations are based on 10 shooters with four shots each. For reference, 
the OPT policy with two shots each is approximately as effective as GREEDY/RANDOM with 
four shots. In this case, optimal use of information is approximately equivalent to doubling the 
firing assets available. If such results survive testing in more realistic scenarios, internetting of 


fires will truly be worth the cost. 


24 





Appendix 


_ This Appendix contains the GAMS file to implement NLP1 with random scenario data. It 
can be amended to handle a given scenario by using input files to read the scenario data. 


$TITLE SHOOTER-TARGET ASSIGNMENT PROBLEM, MINCOST WITH 
EFFECTIVENESS THRESHOLD 


SOFFUPPER ONSYMLIST OFFS YMXREF 


SINLINECOM { } 

option iterlim = 10000000 ; 
option reslim = 10000000 ; 
option minlp = sbb ; 
options optcr = .4; 


SONTEXT 


Original: 5/01/2002 

Authors: Dave Olwell and Alan Washburn 

Revised: yes, 5/09/2002 

Description: Mixed Integer Non-Linear Program designed 
to find optimal assignment of shooter-target 
combinations to minimize cost( be it weight, cube, 
$$, or utils) while reducing the expected fraction 
of remaining target value below a given threshold. 


Uses random assignment of pks, values, and costs. 


Designed to explore characteristics of algorithm prior 
to developing version that imports data from files. 


This is a 100 shooter, 200 target version; 
size can be adjusted by changing indices for shooters and 


target. 
Sponsor: Mr. Mike Bauman, TRAC. 


S$OFFTEXT 


* Option memnodes = 100000; {recommended by SBB documentation, but not allowed by 
GAMS} 
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Sets 
s shooters / sht1*sht100 / 
t targets /tgtl*tgt200/; 
scalar max rounds maximum number of rounds fired by a shooter at a target / 3 / 
b maximum fraction remaining  /.3/ 
max Vv maximum target value /10/ 
max_c maximum shooter cost /8/ 
good q initial search bound 
good _c initial search bound 
good _v initial search bound 
ttiv total initial target value (computed below) 
good _gq=.3; 
good_c=ceil(max_c/2); 
good _v=ceil(max_v/2); 
nnn naan PARAMETERS ---------------------------22--- 22-22-22 nee == 
Parameter 


c(s) cost to fire a round from shts; 

c(s) = floor(uniform(1,max_c+ 1)); 
{ c(s) is uniformly distributed over (1,2,3... max_c)} 
Parameter 

v(t) value of target t ; 

v(t) = floor(uniform (1,max_v + 1)); 
{ v(t) is uniformly distributed over (1,2,3... max_v)} 
Parameter q(s, t) miss probabilities; 


q(s,t) = max(floor(uniform(1,10))/10, floor(uniform(.5,1.5))); 


{ 50% of the qs are 1 by random assignment here, 
rest are uniform over (.1, .2, .3, ..., .9)} 


{ c(s), v(t), b, and q(s,t) will eventually be read in from a data file } 


ttiv = sum(t, v(t)); 
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*__------------ VARIABLES --------------- 
variables 
x(s, t) rounds from shti to tgtk 
Z total cost of firing solution 
ttev expected remaining target value 
frac ttev divided by ttiv 


truecost cost without non-integer penalty; 


integer variable x; 


2 8 A ee 8 2 A ae 2 2 A A OE EAE TH 41.9] bounds 2 2 3 2k 2k 2k 2 2 ofc 2k 2k ik oe 2 2K 2 2K 2c Kk 2K 
x.up(s,t) = max_rounds; 

{will eventually be bounded via input file , 

only fire MAX ROUNDS For a given (s, t)} 


loop((s,t), 1f( q(s,t) = 1, x.up(s,t) = 0)) ; 
{bound the solution space reasonably, can't fire at guaranteed 
miss} 


oe 2 ee fe he ae 2 AK ie 2 2 eK ae 2 EAE EET It 9] variable values* * #*F###EK* 
x.l(s, t)$(q(s, t<=good_q and c(s)<= good _c and v(t)>= good_v) = 1; 
{ pick the high value shots initially to seed search space} 


Equations 
cost define objective function 
threshold __ ttev less than b of ttiv 
ttevdef defines expected remaing target value 
fracdef defines frac; 


cost .. Z =e= sum( (s, t) , c(s) * x(s, t)) ; 

threshold .. ttev =l= b * ttiv ; 

ttevdef .  ttev =e= sum(t, v(t)* prod(s,q(s, t) **(x(s, t)))); 
fracdef = frac =e= ttev/ttiv; 


model probde /all /; 


ZT 











probde.scaleopt = 1; 


probde.nodlim = 50000; 
* probde.optfile = 1; 
probde.tryint = .25; 


KKEKKKEEKK First pass -- solve relaxed problem ******* 
oe ee oe 6 2 2 OK KK KK 

solve probde using rminlp minimizing Z ; 

Oi OK 26 OK 26 2 2 2K OK oe 2K OK 


2K 2 2 3K 26 2K 9k 2 2K 2c 2 2 OK 2K 2 OK OK OK formatted output stuff OK 2 Die 2 fe 2 2 2 2h 24 26 2 2 oe 2 oe I 2K Oe ic 2k 2k 2 ie ak 2K OK Ok Ke 


display x.], x.m ; 
file output1 /outputl.csv/, results /results1.txt/; 


put results; 

put @10, "FORMATTED RESULTS ' ///; 

put ‘Date of run',@15,system.date ; 

put /,'Time of run',@15,system.time; 

put /, card(s), ' shooters and ', card(t), ' targets’ ///; 


put ‘Relaxed cost',@15,z.1/; 

put 'Rounded cost',@15, sum((s, t), c(s)*round(x.l(s, t))) /; 

put "Rounded effectiveness’/; 

put @15, sum(t, v(t)/ttiv * prod(s,q(s, t) **(round(x.|(s, t))))) ///; 
put 'Non-zero shooter-target pairs and round count'//; 


put 'Shooter',@12,'target',@24,'Rounds' /; 

loop((s, t), if(x.1(s, t)>0,put s.t],@12,t.tl@24,x.l(s, t):5:1/) ; 
put //,'cost'’, @10, z.1 :6:1 ///; 

put ‘shooter totals’ /; 

loop(s,put /,s.tl, @10, sum(t,x.1(s, t)):6:1 ); 


put ///, 'ttiv', @25,ttiv:10:4; 

put /, 'ttev', @25,ttev.1:10:4; 

put ///, "B', @25, b:4:3 ; 

put /,'fraction surviving’, @25, frac.1:6:5 ; 


put ///, "List of costs’; 


loop(s, put /, s.tl, @10, c(s)); 
put//, 'List of target values'; 
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loop(t, put /, t.tl,@10,v(t)); 
putclose results; 


put output1l; outputl.pc=5; outputl.pw=10000; 
put ‘target/shooter qs’; 

loop(s, put s.tl); 

loop(t, put / t.tl, loop(s, put q(s, t))); 

put ///; 

put ‘target/shooter rounds’; 

loop(s, put s.tl); 

put /,'totals', loop(s, put sum(t,x.1(s, t))); 
loop(t, put / t.tl, loop(s, put x.1(s, t))); 


putclose output]; 


oF 3 2 2K 2K 2k 


loop((s, t), x.up(s, t)$(x.l(s, t)}=0)=0); 
loop( (s, t), x.l(s, t) = round(x.1(s, t))); 


KEEKEEKLELEES econd pass - solve integer problem with loose optcr ***** 
OF 46 2 26 2K 2c 2K 26 2 ok 2K Ok 


solve probde using minlp minimizing z ; 
6 2h 2c 2c 6 24 3 2c 246 oie ok ok 


Teeter Toren nee es OULDUT SUIT eee Ce On ee ee eee eke ee ee ee 


display x.1, x.m ; 
file output2 /output2.csv/, results2 /results2.txt/; 


put results2; 

put @10, "FORMATTED RESULTS' //; 

put "Date of run',@15,system.date ; 

put /,"Time of run',@15,system.time; 

put /, card(s), ' shooters and ', card(t), ' targets’ /; 
put ‘True integer cost',@15,z.1 //; 


put 'Effectiveness',@20,frac.1:5:4 ///; 
put 'Non-zero shooter-target pairs and round count'//; 


put 'Shooter',@12,'target',@24,’"Rounds' /; 
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loop((s, t), if(x.l(s, t)>0,put s.tl,@12,t.t],@24,x.I(s, t):5:1/) ; 
put //,'cost', @10, z.1 :6:1 ///; 

put 'shooter totals’ /; 

loop(s,put /,s.tl, @10, sum(t,x.1(s, t)):6:1 ); 


put ///, 'ttiv’, @25,ttiv:10:4; 

put /, 'ttev', @25,ttev.1:10:4; 

put ///, "B', @25, b:4:3 ; 

put /,'fraction surviving’, @25, frac.1:6:5 ; 


put ///, 'List of costs’; 

loop(s, put /, s.tl, @10, c(s)); 
put//, 'List of target values'; 
loop(t, put /, t.tl,@10,v(t)); 
putclose results2; 


put output2; output2.pc=5; output2.pw=10000; 
put 'target/shooter qs’; 

loop(s, put s.tl); 

loop(t, put / t.tl, loop(s, put q(s, t))); 

put ///; 

put 'target/shooter rounds’; 

loop(s, put s.tl); 

put /,'totals', loop(s, put sum(t,x.1(s, t))); 
loop(t, put / t.tl, loop(s, put x.l(s, t))); 

putclose output2; 
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